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ABSTRACT 

We model the formation of black hole-neutron star (BH-NS) binaries via dynamical 
interactions in globular clusters. We find that in dense, massive clusters, 16-61% of 
the BH-NS binaries formed by interactions with existing BH binaries will undergo 
mergers driven by the emission of gravitational radiation. If the BHs are retained 
by the cluster after merging with a NS, the BHs acquire subsequent NS companions 
and undergo several mergers. Thus, the merger rate depends critically upon whether 
or not the BH is retained by the cluster after the merger. Results from numerical 
relativity suggest that kick imparted to a ~ 7 Mq BH after it merges with a NS will 
greatly exceed the cluster's escape velocity. In this case, the models suggest that the 
majority of BH-NS mergers in globular clusters occur within 4 Gyrs of the cluster's 
formation and would be unobservable by Advanced LIGO. For more massive BHs, on 
the other hand, the post merger kick is suppressed and the BH is retained. Models 
with 35 Mq BHs predict Advanced LIGO detection rates in the range 0.04 — 0.7 yr" 1 . 
On the pessimistic end of this range, BH-NS mergers resulting from binary-single star 
interactions in globular clusters could account for an interesting fraction of all BH-NS 
mergers. On the optimistic end, this channel may dominate the rate of detectable 
BH-NS mergers. 

Key words: black hole physics - gravitational waves - stars: neutron - globular- 
clusters : general 



1 INTRODUCTION 

Surveys of the Milky Way globular cluster (GC) system 
have discovered 13 low mass X-ray binaries (Deutsch et al. 
2000; Sidoli et al. 2001) and 143 pulsars, many of which are 
recycled, millisecond pulsars (MSPs) and/or members of 
exotic binaries (Lynch & Ransom 2011; Freire et al. 2007; 
Hessels et al. 2006; Ransom et al. 2005; D'Amico et al. 
2001). These observations reveal that 1) some fraction of 
the neutron stars (NSs) formed in GCs are retained de- 
spite their natal kicks, and 2) multibody interactions in 
the cores of GCs greatly enhance the formation rate of bi- 
nary systems that contain a NS, resulting in an over abun- 
dance of low mass X-ray binaries and MSPs relative to the 
field (Katz 1975; Clark et al. 1975; Verbunt & Hut 1987). 
It has long been recognised that binary-single star interac- 
tions between NSs and primordial binaries can result in the 
NS exchanging into the binary by ejecting one member and 
becoming bound to the other (Hills 1976; Hut & Bahcall 
1983; Sigurdsson & Phinney 1995; Ivanova et al. 2008). 
Such binary-single star interactions could also result in the 
formation of black hole (BH)-NS binaries (Sigurdsson 2003; 
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Devecchi et al. 2007), which are expected to produce gravi- 
tational waves detectable by the ground-based interferome- 
ter LIGO if they coalesce (Abbott et al. 2009). Observations 
of a BH-NS merger would provide a test of General Rela- 
tivity in the dynamical strong-field regime and could probe 
the NS equation of state (Kyutoku et al. 2009; Duez et al. 

2010) . Additionally, BH-NS mergers are also a likely pro- 
genitor of some short gamma-ray bursts (Nakar 2007). 

Whether or not dynamical interactions in GCs can ef- 
ficiently produce BH-NS binaries is unclear. Models sug- 
gest that the ~ 10 — 300 BHs that form within a GC 
will rapidly sink to the core of the cluster and inter- 
act with, and eject one another, severely depleting the 
cluster's BH population within 1 Gyr (Kulkarni et al. 
1993; Sigurdsson & Hernquist 1993; O'Leary et al. 2006; 
Moody & Sigurdsson 2009; Banerjee et al. 2010; Aarseth 
2012). Even though most of GCs BHs are ejected early in its 
evolution, there are models (Miller & Hamilton 2002) and 
mounting observational evidence that show some GCs retain 
at least one BH. Several extragalactic GCs harbour highly 
variable X-ray sources with luminosities > 10 39 erg s _1 , well 
in excess of the Eddington luminosity of an accreting NS, 
that are strong candidates for BHs in GCs (Maccarone et al. 

2011) . Even if a GC only retains a single BH, it may be 
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possible for this BH to merge with multiple NSs if the 
BH exchanges into another NS-containing binary after each 
merger. 

What is the maximum rate for such successive merg- 
ers? If we assume that the BH and NS merge instanta- 
neously, then the maximum merger rate is equal to the rate 
at which the BH can exchange into a NS-containing binary, 
Rex = ^NSCexWrei) where tins is the density of NS-binaries, 
(j cx is the cross section for an exchange, and v Te i is the rela- 
tive velocity of the BH and the binary. Ivanova et al. (2008) 
found that GCs tend to retain ~ 220 NS per 2 x 1O 5 M . 
Due to mass segregation, NSs will be concentrated near the 
centre of a GC, and the central NS fraction, /ns, can be as 
high as 0.1. Assuming that all of the NSs at the centre of the 
cluster are in a binary and that the collision cross section 
is dominated by gravitational focusing, we can approximate 
the maximum BH-NS merger rate as 
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Here, n is the total central stellar density, a is the typical 
semi-major axis of a NS binary, Mbh is the BH mass, and 
v m is the mean central velocity dispersion. 

In this paper we present detailed models that investi- 
gate the rate at which the BHs and NSs retained by a GC 
will form binaries that merge via the emission of gravita- 
tional radiation and the corresponding detection rate for 
LIGO. We will begin by discussing the methods used in our 
simulations in Section 2. In Section 3.1, we will describe the 
results of simulations that only allowed BHs to merge with 
a single NS, and in Section 3.2 we will describe how the re- 
sults change when the BH is allowed to under go successive 
mergers as sketched above. We will conclude in Section 4 
by computing the LIGO detection rate and comparing our 
work with previous calculations. 



2 METHOD 

We modelled the formation of BH-NS binaries through 
binary-single star interactions using the method presented 
in Sigurdsson & Phinney (1995) and Mapelli et al. (2005). 
For each simulation, an ensemble of binaries was evolved in 
a static background cluster. The background clusters were 
multi-mass King models, with the stars of each mass group 
a following the isotropic distribution function 



o 



1] e > 
e < 
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where v ma is the mass group's core velocity dispersion, no Q 
is a normalising constant, and e = $>(r) — vt/2 is the relative 
energy per unit mass. Here, v* is the velocity and ^(r) is 
the relative potential given by $(r t ) — $(r), where $(r) is 
the gravitational potential with respect to infinity at a dis- 
tance r from the cluster centre and rt is the cluster's tidal 
radius. To explore how the BH-NS merger rate is impacted 
by the structural properties of the GC, we varied the cen- 
tral stellar number density, n c , the number density weighted 
mean central velocity dispersion, v m , and the King parame- 
ter Wo = $f(0)/v m to generate clusters of different mass and 



create a set of background clusters with core radii r c and 
concentrations c = \og(rt/r c ) similar to observed clusters. 
Furthermore, these clusters probe a range of core interaction 
rates V oc 7ic' 5 r;?, a quantity that is known to correlate with 
another product of dynamical interactions, the number of 
X-ray binaries present in a cluster (Pooley et al. 2003). We 
list the values used in our models in Table 1. 

For each cluster, we used the initial stellar mass function 



£(m) oc 



m < 0.35 M 
m > 0.35 M 



(3) 



(Kroupa 2001), and binned the evolved population into 10 
mass groups. The main sequence (MS) turnoff mass, m% , 
was set to 0.85 Mq and stars with initial mass above the 
turnoff mass were assumed to be completely evolved. The 
evolved stars fell into one of three groups, those with initial 
mass rrii < 8 Mq evolved into white dwarfs (WDs) with 
masses given by mwD = 0.38 + 0.12 rm (Catalan et al. 
2008). Stars with initial mass in the range 8 Mq < rm < 
25 Mq were assumed to form 1.4 Mq NSs and those with 
initial mass > 25 Mq formed BHs. 

We performed calculations using two different values 
for the BH mass. Recent statistical analyses of low mass 
X-ray binaries in the Galaxy have found that the masses 
of BHs are narrowly distributed around 7 Mq (Ozel et al. 
2010; Farr et al. 2011). Motivated by this, we have run mod- 
els in which all BHs have a mass of 7 Mq. However, the 
progenitors of the field BHs considered in the above studies 
were stars with much higher metallicity than those found in 
GCs. Lower metallicity stars have reduced mass loss rates 
and could evolve to BHs with masses as high as 80 Mq 
(Fryer & Kalogera 2001; Fryer et al. 2002; Belczynski et al. 
2010). There is also observational evidence that BHs in GCs 
are more massive than those in the Galaxy. Properties of 
the ultraluminous X-ray sources observed in GCs associ- 
ated with NGC 4472, M31, and NGC 1399 are consistent 
the presence of a BH of > 30 Mq (Maccarone et al. 2007; 
Barnard et al. 2011; Irwin et al. 2010; Clausen et al. 2012). 
Accordingly, we also ran models with 35 Mq BHs to allow 
for the possibility that GC BHs are more massive than those 
in the field. 

We assumed that the clusters retained 20% of the NSs 
that were formed within them. We further assumed that 
nearly all BHs formed in the cluster were promptly ejected 
by self interaction, leaving 0, 1, or 2 BHs (Kulkarni et al. 
1993; Sigurdsson & Hernquist 1993; O'Leary et al. 2006; 
Moody & Sigurdsson 2009; Banerjee et al. 2010; Aarseth 
2012). We mimic this by truncating the high mass end of 
the initial stellar mass range so there is a single BH in the 
GC model. We are, therefore, modelling the late evolution 
of the cluster, and hence recent times, so the merger rates 
calculated here apply to the local universe. In all of our 
simulations, we assumed that the remaining BH was in a 
binary. The properties of such a binary will have undergone 
extensive modification driven by both stellar and dynami- 
cal evolution. Consequentially, the initial configurations as- 
signed to the BH binaries in our models represent the result 
of complicated processes that we do not attempt model. 
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Table 1. Background Globular Cluster Models Parameters 



Model 


n c 


v m 


W 


M 


c 


Name 


(pc- 3 ) 


(km s- 1 ) 




(Mq) 




A 


1 x 10 4 


6 


6 


1.0 x 10 s 


1.20 


13 


1 x 10 5 


10 


10 


5.2 x 10 5 


1.71 


C 


5 x 10 5 


11 


13 


7.2 x 10 5 


1.93 


D 


1 x 10 6 


12 


15 


1.1 x 10 6 


2.06 



2.1 Dynamics 

For each simulation, we evolved 2000 binaries, one at a time, 
in one of the background clusters described above. Each 
binary contained a BH and a companion drawn at ran- 
dom from the cluster's evolved mass distribution. The initial 
semi-major axis, a, was drawn from a distribution that is fiat 
in log a between 10 -3 au and a max (Abt 1983). We used a 
different value of a max for each background cluster, namely 
100 au, 33 au, 15 au, and 10 au for clusters A, B, C, and 
D, respectively. These values of a max were chosen to cover 
the range of possibilities that might result from the compli- 
cated evolutionary histories of the BH binaries mentioned 
above. The initial eccentricity was drawn from /(e) oc 2e 
(Duquennoy & Mayor 1991). 

Each binary was initially placed in the cluster core, with 
its position selected from the radial density distribution of 
the third most massive mass group. The results are insen- 
sitive to this somewhat arbitrary choice because the initial 
position and velocity of the binary are rapidly forgotten. 
The binaries were evolved in the cluster potential with dy- 
namical friction and random kicks to account for scattering 
by individual stars. At each point along the binary's trajec- 
tory, we computed the probability that it would experience 
a strong encounter with a single star using 



= At ^ J n c 



(r)a(v, v*)\v — v f \f a (v*)d 3 v* , (4) 



where the subscripts a correspond to the mass groups, At 
is the time step, n is the number density, f(v t ) is the ve- 
locity distribution of single stars given in equation 2, v and 
v, are the velocity of the binary and a single star, respec- 
tively, and a is the encounter cross section. Throughout this 
paper we are interested in encounters that alter the energy, 
angular momentum or components of the binary. Previous 
studies have shown that the bulk of these changes involve 
a background star approaching the binary centre of mass 
closer than a few times the binary's physical size and con- 
versely that numerous, more distant encounters have rel- 
atively small effect. We choose the critical pericentre dis- 
tance for simulated encounters to be p = a[4 + 0.6(1 + e)] 
where a is the binary semi-major axis and e is its eccentric- 
ity (Hut & Bahcall 1983). The corresponding cross section 
for such encounters is 



, , 2 2irG(m bin + m a )p 
cr{v,v*) = np H ■ 12 , 



where mbm is the mass of the binary. An encounter took 
place if Penc was greater than a random number drawn 
from a uniform distribution between and 1. Given that 
an encounter occurred, we randomly chose which particular 
mass group was involved based on the fractional encounter 



probabilities of all the groups. Then, the three-body interac- 
tion between the stars was integrated explicitly as described 
in Sigurdsson & Phinney (1993). There were many possible 
outcomes to the three-body encounter, including mergers, 
exchanges, and disruption of the binary, but if a binary still 
existed at the end of the interaction, then we continued to 
evolve the binary in the cluster. 

We performed two types of simulations. In the first class 
of simulations we were only concerned with the fate of the 
initial BH binary, while in the second class we sought to 
determine the ultimate fate of the BH. In the former set, a 
run continued until either the initial binary was disrupted, 
merged, or ejected, or the run had covered 10 10 yrs. In the 
latter set of models, we continued a run when the merger 
or disruption of a binary resulted in an isolated BH (see 
Section 2.3). These runs ended only after the BH was ejected 
from the cluster or the maximum run time of 10 10 yrs was 
reached. In both simulation types, a small fraction of the 
runs ( < 5%) were aborted because integration of a three- 
body interaction exceeded the maximum number of allowed 
steps, 8 x 10 7 . 



2.2 Gravitational radiation effects on the 
evolution BH— NS binaries 

Between interactions the eccentricity and semi-major axis of 
each BH-NS system were decreased to account to account 
for the emission of gravitational radiation using the rela- 
tions given in Peters (1964). For the simulations in which 
we used a BH mass Mbh = 7 Mq, we stopped the calcula- 
tion when a BH-NS binary coalesced. Numerical relativity 
simulations have found that the remnant formed by a BH-NS 
merger at this mass ratio will receive a kick from anisotropic 
emission of gravitational radiation that exceeds the GC's 
escape velocity (Shibata et al. 2009; Etienne et al. 2009; 
Foucart et al. 2011). In the models where we used Mbh = 
35 Mq, the post-merger recoil velocity was suppressed be- 
cause of the much smaller mass ratio (q — M^s/Mbh)- 
For the case of non-spinning point masses, the merger of 
a 1.4 Mq NS and a 35 Mq BH would result in a recoil ve- 
locity of 14.6 km s , well below the escape velocity for the 
clusters models that we considered (Gonzalez et al. 2007). 
However, we allowed for spinning BHs and NSs, and com- 
puted the recoil velocity using the parameterisation given 
in Campanelli et al. (2007, however, see Hirata 2011 for dis- 
cussion of a resonant recoil, not included in these parame- 
terisation, that can dominate when the BH spin is nearly 
maximal and q <C 1). At the beginning of each run, the BH 
spin parameter was randomly selected from a uniform distri- 
bution between and 1. When a BH-NS binary merged, the 
NS spin was randomly selected from a uniform distribution 
between and 0.7 (a spin parameter of 0.7 corresponds to 
a NS spinning at the break up frequency Miller et al. 2011) 
and the two spins were randomly oriented with respect to 
one another and to the orbital plane. The result was a range 
of kick velocities (~ 10 — 100 km s _1 ), some of which were 
low enough that the single BH was retained by the cluster 
after the merger. 
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2.3 Formation of a new binary from a single BH 
by stellar interactions 

In some of our models, when the BH lost its binary compan- 
ion, either through an interaction with a background star or 
as the result of a gravitational radiation driven merger with 
a NS, we continued to evolve the now single BH in the clus- 
ter. As the BH moved through the cluster, we sought to 
determine the single BH's next binary companion by gen- 
crating binaries with a semi-major axis a, an eccentricity e, 
a primary of mass mi, and a secondary of mass 7712. The 
semi-major axis and eccentricity were drawn from the same 
distributions as those of the initial BH binaries. The only 
difference is that we have chosen a max = 1 au for these sys- 
tems to ensure that the binaries are hard. Our choices for 
mi and 7712 are described below. 

We determined whether or not the single BH would 
interact with the binary by temporarily replacing the sec- 
ondary with the BH. We then calculated the probability 
that this temporary binary would interact with a back- 
ground star of mass 7712 (the secondary), -P(i,bh)+2, us- 
ing equation 4. Next, P(i i bh)+2 was scaled to -P(i,2)+bh by 
comparing the rates of two related exchange encounters, 
[(l,BH) + 2-> (2,BH) + 1] and [(l,2)+BH-> (l,BH) + 2]. 
The rate for an arbitrary interaction between a binary (a,b) 
and single star c can be expressed in the form 

. , Gm a m b (m a + m b + m c ) a in na(X) 

K y x > = ~i — m ' w 

m c (m a +mb) v 

where X refers to the type of interaction (e.g., one of the 
desired exchange reactions), v is the relative velocity of the 
binary and the single star, is the semi-major axis of 
the binary, n is the density of single stars, cr(X) is a nu- 
merically determined, dimensionless cross section, G is the 
gravitational constant, and rru is the mass of star i. Using 
equation 6, we computed the relative rates of the above ex- 
change interactions, assuming that ai n , n, and v were the 
same in each case. 

The required dimensionless cross sections were 
computed by Sigurdsson & Phinney (1993), who found 
er[(l,BH) + 2 -> (2,BH) + 1] ~ 2 x er[(l,2) + BH -> 
(1,BH) + 2]. In addition, we made the simplifying approxi- 
mation that mi = 7712 <S Mbh, because the BH was much 
more massive than either the primary or the secondary. This 
reduced the ratio of the mass factors to 7712 /2A/bh . Combin- 
ing these, we approximated the probability that the single 
BH would interact with a binary containing the primary and 
the secondary as 

P h. 7712 P (7\ 

Ml,2 +BH ~ z J-Y~r r (l,BH)+2. {<) 

1 — jb Mbh 

where fb is the binary fraction and relates the densities of 
single and binary stars. Note that the right hand side of 
this expression diverges as fb — > 1. Thus, for large values of 
fb, use of equation 7 will overestimate the rate at which a 
single BH can exchange back into a binary. Furthermore, 
we have used the relative rates of two exchange encoun- 
ters to scale the total encounter probability. The exchange 
cross section is a significant fraction of the total interaction 
cross section, but we will explore the consequences of this 
choice in Section 4.1. We chose to compute the encounter 
probability using the scaling relation given in equation 7 
because the distribution of binaries in globular clusters is 



poorly understood and this technique offers an efficient and 
adequate estimation of the rates within the limitations of 
our code. When an encounter was deemed to have occurred, 
we restored the original binary and explicitly integrated the 
three-body interaction between the binary and the BH, ac- 
cepting any outcome. We repeated this process until the BH 
exchanged into a binary, the BH was ejected from the clus- 
ter, or the maximum integration time was reached. If the 
BH exchanged into a binary, we continued the simulation as 
described in Section 2.1. 

We explored several different binary populations and 
values of fb. In one set of models, we drew the primary's 
mass from the cluster's evolved mass distribution and re- 
quired that mi > m to , but we did not chose the secondary's 
mass, initially. Instead, we used a binary consisting of the 
primary and the BH to compute the scaled probability that 
a BH would interact with binary containing the primary 
and a member of any mass group. If an encounter occurred, 
then the secondary was chosen to belong to the mass group 
that interacted with the temporary binary. We then inte- 
grated the three-body interaction as described above. This 
method of selecting the secondary resulted in a binary pop- 
ulation that is comprised largely of systems that contain a 
NS. For the simulations in cluster D, more than 50% of the 
binaries that interacted with a single BH contained a NS. 
We have labeled the runs that use this binary population 
OPT, for optimised, because this method produces a binary 
population that is tuned to interact with the single BH. Fur- 
thermore, because many of the binaries in this population 
already contain a NS, the single BH can gain a NS compan- 
ion immediately. 

The second binary population that we considered was 
closely linked to observational constraints. A number of 
studies have used the fact unresolved MS binaries appear 
brighter and redder than single MS stars to investigate bi- 
naries in GCs. The binary fraction can be measured by 
comparing the number of objects in the offset binary se- 
quence to the number of stars that lie along the MS, 
with observed values typically found in the range f b = 
0.05 — 0.3 (e.g., Romani & Weinberg 1991; Sollima et al. 
2007; Davis et al. 2008; Milone et al. 2012). We generated 
a population that resembles these observed populations by 
independently drawing two stars from the cluster's evolved 
stellar mass distribution and requiring that one of the stars 
have a mass m > m to . The result was a collection of sys- 
tems that consisted primarily ( > 85%) of MS-MS binaries. 
We have labeled simulations that use this binary population 
OBS, and considered values for the binary fraction within 
the observed range: 0.05, 0.1, and 0.2. 

We also considered a population that was motived by 
detailed models of binary evolution in GCs. Analytic, Monte 
Carlo, and N-body models have indicated that over time 
the hard binary fraction will increase in the cluster core 
(Sollima 2008; Fregeau et al. 2009; Hurley et al. 2007). For 
such models to be consistent with the observations described 
above, the initial binary fraction in GCs would have to have 
been quite low. Alternatively, as Fregeau et al. (2009) sug- 
gested, if many of the binaries in a GC consisted of a MS 
star and a WD or two WDs, they would be "hidden" from 
observers. Following on this result, we have run simulations 
that use a binary population similar to that described in 
table 1 of Fregeau et al. (2009). This population consisted 
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Figure 1. Initial locations of the BH-NS binaries in the (1 — e 2 )- 
a plane. The points shown here are from simulations in clus- 
ter model D with A/bh = 7 Mq. tGW is constant along the 
solid lines. From left to right, the lines correspond to tGW = 
10 5 , 10 6 , 10 7 , 10 8 , 10 9 , 10 10 yrs. 



of 44% MS- WD binaries, 32% WD- WD binaries, 23% MS- 
MS binaries, and 1% binaries containing a NS. Runs that 
used this binary population were labeled FIR. Furthermore, 
because only 23% of this binary population would be mea- 
sured by the observations described above, we have used 
larger values of with this population: 0.2, 0.5, and 0.75. 



3 RESULTS 

In all of our simulations, binary-single star interactions pro- 
duced BH-NS binaries. However, the number of BH-NS 
formed and the properties of these systems were very sen- 
sitive to the structure of the GC and the assumed binary 
population. Our simulations are summarised in Table 2. The 
initial conditions for each simulation, including the cluster 
model, the binary population (where applicable), the binary 
fraction (where applicable), and Mbh, are given. Also listed 
are /bhns, the fraction of runs that produced a BH-NS bi- 
nary, JVbhns, the total number of BH-NS binaries formed 
in all 2000 runs, -/V mrg , the number of BH-NS mergers dur- 
ing the simulation, t cx , the average amount of time taken 
for a single BH to exchange back into a binary, iV max , the 
maximum number of mergers in a single run, and IZgc, the 
BH-NS merger rate for that model. We describe the results 
in detail below, beginning with models that assumed that 
the BH was ejected from the GC after merging with a NS, 
followed by descriptions of models with more massive BHs, 
which allowed for a single BH to merge with multiple NS. 



Figure 2. The actual amount of time it took a BH-NS binary 
to merge £qw vs - the amount of time expected for the BH-NS 
binary to merge, given the binary's initial configuration, 4gw 
for the BH— NS binaries that merge in a simulation of cluster 
model D with Mbh = 7 Mq. Note that many systems formed 
with tGW a Hubble time. Subsequent interactions between 
these binaries and single stars greatly reduced (qw We excluded 
systems where the initial binary was randomly chosen to contain 
a BH and NS in our rate estimates to avoid initial condition 
bias and because these would happen at high rcdshift and be 
unobscrvable. 



3.1 Single merger scenario 

3.1.1 Without single BHs 

In this section we describe simulations that ended if the 
binary containing the BH was disrupted or merged. These 
models did not require any assumptions about the GCs bi- 
nary population as a whole because they were terminated 
as soon as BH lost its companion. We have computed one 
such simulation in each of our background clusters, and in 
all of these simulations we assumed Mbh = 7 Mq. 

In the lowest density cluster considered here, cluster A, 
only eight BH-NS binaries were formed in our 2000 runs, 
none of which merged during the 10 10 yr simulations. In 
all of the other clusters, however, some fraction of the BH- 
NS binaries merged during the simulation. Figure 1 shows 
the initial semi-major axis and eccentricity of the BH-NS 
systems formed in simulations using cluster model D. Also 
plotted are lines of constant gravitational radiation merger 
time. A majority of the systems lie above the right most line, 
which corresponds to a merger time 10 10 yr. If these systems 
were to evolve in isolation, it would take more than a Hubble 
time for them to merge. However, in the massive GCs consid- 
ered here subsequent interactions with single stars extracted 
energy from the binary's orbit and drove some of these sys- 
tems to high eccentricity and/or small separations, greatly 
accelerating the merger process. This mechanism is illus- 
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trated in Figure 2. BH-NS binaries with expected merger 
times, given their initial eccentricity and semi-major axis, 
tewo > 10 7 yr underwent several interactions with single 
stars after their formation, which altered the orbital pa- 
rameters of the binary and reduced the actual amount of 
time between the formation of the BH-NS binary and the 
merger of the binary, £gw, by several orders of magnitude. 
In these simulations, binary-single star encounters not only 
produced BH-NS binaries, but in some cases enhanced the 
rate at which these systems merged. 

In other situations, subsequent interactions reduced the 
merger rate. In models that used cluster B, many BH-NS 
binaries were destroyed when a WD exchanged into the bi- 
nary. In this relatively low concentration cluster, the density 
of WDs was much greater than the density of NSs, so mul- 
tiple encounters with WDs overcame the low probability of 
a WD exchanging into the BH-NS binary. This property of 
low concentration clusters has implications for the forma- 
tion of mass transfer binaries, and hence X-ray luminosity 
of such clusters, which we will investigate in a future publi- 
cation. At the end of most runs in cluster B, the BH remains 
in the cluster. The BH is ejected in 3% of the runs. Of those 
ejected, 42% are ejected after a BH-NS merger and 58% 
are ejected after a super-elastic collision between the binary 
and a background star. In the majority of runs (63%), the 
BH remains in the cluster with a non-NS binary companion. 
The BH loses its companion but remains in the cluster in 
the remaining 34%. 

This over-density of WDs did not occur in the high 
concentration clusters, C and D. In these simulations, if a 
BH-NS binary was destroyed by an exchange, it was usually 
another NS exchanging into the system. The simulations 
suggest that in the clusters with n > 5 x 10 5 pc~ 3 , there is 
a high likelihood that a BH retained by the cluster after the 
initial phase of self interaction will become part of a BH-NS 
binary that merges within a Hubble time. This, coupled with 
the high post-merger kicks predicted by numerical relativity, 
suggests that if GC BHs have the same mass distribution 
as those in the field, BHs will only be retained in massive 
clusters with n < 10 5 pc~ 3 . 

The distribution of BH-NS merger times, tmrg, for 
the different GCs models are shown as solid histograms in 
Figure 3. We note that t m rg > tew because t mrg includes the 
time taken to produce the BH-NS binary. Encounters oc- 
curred frequently in simulations of the more massive, higher 
density clusters, so the BH-NS binaries were formed earlier 
and, as described above, encounters after their formation ac- 
celerated the merger process. For each cluster we calculated 
the merger rate as TLgc = N m rg/(Nt m ax), where N is the 
total number of binaries simulated and t m ax = 10 Gyr is 
the length of each run. In each cluster, N was slightly less 
than 2000 because we discarded the runs where the initial 
binary was randomly chosen to contain a BH and NS. Al- 
though these simulations predicted many BH-NS mergers, 
a substantial faction of the runs ended prematurely because 
the BH lost its companion. For the runs in clusters C and 
D, nearly 40% of simulations ended for this reason. In the 
next section we will examine how the merger rate changes if 
such single BHs are allowed to interact with and exchange 
back into a binary. 



3.1.2 With single BHs 

Here, we describe a set of simulations that allowed BHs that 
lost their companions as the result of a binary-single inter- 
action to continue their evolution in the cluster. Again, we 
have run one such simulation in each of our background clus- 
ter models and assumed that Mbh = 7 Mq . We note that in 
these simulations, each BH could only undergo one merger 
with a NS because with Mbh = 7 Mq the post merger 
kick imparted to the BH exceeded the escape velocity of 
the cluster. In all four simulations, we used the FIR binary 
population with ft = 0.5. 

These changes had little effect on simulations of clus- 
ter A. The interaction rate in cluster A is lower than the 
other clusters, so the initial BH binaries undergo far fewer 
encounters. Only 57 of 2000 runs in this simulation pro- 
duced a single BH, and none of these single BHs exchanged 
back into a binary. On the other hand, allowing the single 
BHs to find new companions increased the number of BH- 
NS binaries that formed and the number of these systems 
that merged in clusters B, C, and D. The average amount 
of time that it took for a BH without a companion to ex- 
change into a binary, t ox , was 3.64, 3.31, and 1.88 Gyrs for 
clusters B, C, and D, respectively. In most cases, a single 
encounter did not result in the BH exchanging into the bi- 
nary, so t cx is longer than the interaction timescale. When 
the BH did exchange into the binary, the resulting systems 
(in all three cluster models) had an average semi-major axis 
of 2.8 au and average eccentricity of 0.8. Because the inter- 
action cross section for such wide and eccentric binaries is 
large, the newly formed systems quickly interacted with ad- 
ditional background stars. In the simulations of cluster D, a 
BH-NS binary was formed in 87% of the runs, and 93% of 
the runs that produced such a binary ended with a BH-NS 
merger. A single run could produce multiple BH-NS bina- 
ries as the result of numerous exchanges, but only a single 
BH-NS merger. For example, in cluster D, the average num- 
ber of BH-NS binaries produced in a single run was 2.1, but 
in one run the BH had ten different NS companions before 
merging. The average BH-NS merger rate in cluster B nearly 
doubled when we allowed single BHs to gain new compan- 
ions. In clusters C and D the value more than tripled. 

The dotted histograms in Figure 3 show the distribution 
of BH-NS merger times in this set of models. These dotted 
histograms are the combination of two t mrg distributions. 
One distribution, similar to the solid histograms, is due to 
BHs that merge with a NS before losing their binary com- 
panion. The second distribution is for mergers that occurred 
after the BH lost its binary companion and had to exchange 
back into a binary before finally gaining a NS companion 
and, eventually, merging. For many of these mergers, t mrg is 
much larger than in the previously discussed cases in which 
the BH merged before losing its original companion. Merg- 
ers that occurred after the BH became single dominate the 
t mrg distributions. So, in addition to an increased average 
BH-NS merger rate, these runs also predict that the bulk 
of these mergers will occur at later times. However, many 
of the mergers, especially those in cluster D, still occurred 
within 3 Gyr of the start of the simulations. This suggests, 
again, that if GC BHs are of similar mass to those found 
in the Galaxy and the large post-merger kicks predicted by 
numerical relativity are correct, most of the BHs that sur- 
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Table 2. Globular Cluster Models and Merger Rates 



Cluster 


Binary 


ft 




/bhns 


Abhns 


N mrg 




-^Y max 




Model 


Population 




(M ) 








(yr) 




(yr- 1 ) 


Single merger scenario without single black holes (Section 3.1.1) 


A 






7 


4 x 10~ 3 


8 












13 






7 


0.080 


202 


28 






1.4 x 10" 12 


C 






7 


0.22 


732 


282 






1.4 x lO" 11 


D 






7 


0.34 


1423 


465 






2.4 x 10~ n 






Single 


merger 


scenario with 


single black holes (Section 3.1.2) 






A 


FIR 


0.5 


7 


5.0 X 10~ 3 


10 





> 10 10 







13 


FIR 


0.5 


7 


0.18 


470 


52 


3.64 x 10 9 




2.7 x 10~ 12 


C 


FIR 


0.5 


7 


0.66 


2568 


1002 


3.31 x 10 9 




5.1 x 10~ n 


1) 


FIR 


0.5 


7 


0.87 


4296 


1626 


1.88 x 10 9 




8.2 x 10~ n 


Multiple merger scenario (Section 3.2) 


A 


FIR 


0.2 


35 


3.0 x 10~ 3 


6 





> 10 10 










FIR 




JU 


14 


ozu 


102 


U.Od X 1U 


i 


O.O X 1U 


C 


FIR 


0.2 


35 


0.42 


1217 


806 


3.98 x 10 9 


2 


4.1 x 10~ n 


1) 


FIR 


0.2 


35 


0.64 


2698 


1635 


3.58 x 10 9 


5 


8.3 x 10~ n 


A 


FIR 


0.5 


35 


2.5 x 10" 3 


5 





> 10 10 








[ > 

]_> 


r lix 


U.u 




U.lO 


ODU 


1 9/1 


q i A v 1 n9 

0. 14 X 1U 


1 


(\ A v 1 C\ — 12 
0.4 X 1U 


C 


FIR 


0.5 


35 


0.54 


1882 


1185 


3.16 x 10 9 


3 


6.1 x 10" 11 


D 


FIR 


0.5 


35 


0.84 


5679 


3311 


2.01 x 10 9 


7 


1.7 x 10~ 10 


A 


FIR 


0.75 


35 


1.5 x 10" 3 


3 





> 10 10 








13 


FIR 


0.75 


35 


0.15 


377 


137 


4.22 X 10 9 


1 


7.0 x 10~ 12 


C 


FIR 


0.75 


.35 


0.68 


3072 


1881 


1.89 x 10 9 


1 


9.8 x lO" 11 


1) 


FIR 


0.75 


35 


0.85 


8664 


5266 


8.13 x 10 s 


9 


2.7 x 10~ 10 


A 


OPT 


0.2 


.35 


3.5 x 10" 3 


7 


1 


> 10 10 


1 


5.0 x 10~ 14 


B 


OPT 


0.2 


.35 


0.16 


364 


146 


4.54 x 10 9 


2 


7.6 x 10~ 12 


C 


OPT 


0.2 


35 


0.57 


2076 


1275 


3.57 x 10 9 


3 


6.5 x lO" 11 


1) 


OPT 


0.2 


35 


0.85 


8258 


4128 


1.84 x 10 9 


9 


2.1 x 10~ 10 


A 


OBS 


0.05 


.35 


2.0 x 10~ 3 


4 





> 10 10 








13 


OBS 


0.05 


35 


0.15 


341 


113 


4.79 x 10 9 


1 


7.3 x 10" 12 


C 


OBS 


0.05 


35 


0.41 


1151 


737 


3.94 x 10 9 


2 


3.7 x 10" 11 


D 


OBS 


0.05 


.35 


0.54 


1760 


1089 


4.42 x 10 9 


3 


5.5 x 10" 11 


A 


OBS 


0.1 


35 


5.5 x 10~ 3 


11 


1 


> 10 10 


1 


5.0 x 10" 14 


13 


OBS 


0.1 


35 


0.13 


319 


133 


4.58 x 10 9 


1 


6.8 x 10~ 12 


C 


OBS 


0.1 


.35 


0.40 


1089 


725 


4.05 x 10 9 


2 


3.7 x 10" 11 


D 


OBS 


0.1 


.35 


0.58 


2054 


1250 


4.17 x 10 9 


3 


6.3 x 10" 11 


A 


OBS 


0.2 


35 


3.0 x 10" 3 


6 





> 10 10 








13 


OBS 


0.2 


.35 


0.14 


337 


113 


4.50 x 10 9 


1 


5.8 x 10~ 12 


C 


OBS 


0.2 


.35 


0.43 


1313 


818 


3.81 x 10 9 


3 


4.2 x 10" 11 


1) 


OBS 


0.2 


35 


0.67 


2980 


1701 


3.56 x 10 9 


1 


8.6 x 10- 11 



Properties of the GC models are given in Table 1, the binary populations are described in Section 2.3, is 
the binary fraction, Mbh is the mass of the BH, /bhns i s the fraction of runs in which a BH-NS binary 
was produced, ./Vbhns i a th e total number of BH-NS binaries formed in the 2000 run simulation, 7V mrg is 
the number of BH-NS mergers during the simulation, t GX is the mean amount of time it takes for a single 
BH to exchange back into a binary, JV max is the maximum number of BH-NS mergers in a single run, and 
TZgc is the BH-NS merger rate. 



vived the initial phase of self-interaction will be ejected from 
the cluster after merging with a NS well before the current 
epoch. 



3.2 Multiple merger scenario 

So far, we have considered scenarios in which the BH re- 
ceives a post-merger kick that exceeds the escape velocity 
of the GC. Under this constraint, each BH can undergo at 
most one merger. If GC BHs are more massive than those 
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Figure 3. Distribution of the merger times for BH-NS binaries in clusters D (left panel), C (centre panel), and B (right panel) with 
Mbh = 7 Mq . These merger times include the time taken to form the BH-NS system through binary— single star interactions. The solid 
histogram is based on runs that extended only until the initial BH binary either merged or was disrupted; it shows the distribution of 
tmrg when a merger occurred. The dotted histograms show the merger time distribution for runs that allowed BHs that lost their binary 
companions to exchange into a new binary. In these runs we have used the FIR binary population with = 0.5. With this change, a 
larger number runs produced a merging BH-NS binary, and the distribution shifts towards longer merger times. The larger merger times 
are due to the significant amount of time it takes for the single BH to exchange back into a binary. In all of the runs shown here, the 
BH is ejected from the cluster after the merger. 



found in field, then it is possible for the BH to remain bound 
to the cluster after the merger, regain a NS companion, and 
merge again. In this section we will discuss models that as- 
sumed Mbh = 35 Mq . The magnitude of the recoil velocity 
is approximately proportional to q 2 , and with a BH mass 
of 35 Mq the recoil velocity becomes comparable to the 
GC escape velocity. Thus, we explored the case in which 
the suppressed kick velocity helped but did not ensure that 
the post-merger BH would be retained by the cluster. The 
fraction of post-merger BHs retained in simulations using 
each background cluster depended on that cluster's escape 
velocity v esc « V2Wv m - Among all the models of the bi- 
nary population and varying binary fraction, in simulations 
using clusters A, B, C, and D, 50%, 23%, 20%, and 10% of 
the BHs that underwent a merger were ejected, respectively. 
Note that only two BH-NS mergers occurred in cluster A, 
so the ejection fraction derived from the models is smaller 
than expected. However, in these simulations many of the 
BHs remained bound to the cluster after merging with a NS. 
A fraction of these single BHs then gained new NS compan- 
ions and underwent multiple mergers. 

The number of times a BH can merge with a NS de- 
pended on how quickly the single BH could exchange into 
a binary after each merger. While the post merger kick was 
not always large enough to remove the BH from the cluster, 
the BH was expelled from the core after most mergers. Be- 
fore it was able to exchange into a binary, it needed to sink 
back into the core where the interaction rates were high- 
est. This process occurred rapidly in all of the simulations. 
The mass segregation timescale is proportional to m/A/sH, 
where m is the average stellar mass in the cluster. The BHs 
in these models were 35 Mq, which gave m/MeH between 
0.014-0.027. Typically, a kicked BH would return to the core 
on a timescale O(10 7 ) yrs. As can be seen from Table 2, this 
timescale is short compared to the total time it takes the 



single BH to exchange into a binary. There was, therefore, 
little difference between the t cyi distributions for BHs that 
became single because of a BH-NS merger and those that 
lost their companions after a binary-single interaction. We 
make no distinction between these distributions in the dis- 
cussion below. 

Figure 4 shows how the distribution of t cx is affected 
by the mass of the BH, the structural parameters of the 
GC, and the binary fraction. All of the simulations shown 
here used the FIR binary population. Changing the mass 
of the BH had little impact on both the distribution of t ex 
and its average value. As expected, i ex decreases with in- 
creasing cluster concentration. Increasing the binary frac- 
tion produced slightly narrower i ox distributions with peaks 
shifted to shorter times. The t ox distributions for simula- 
tions using the OBS binary population with = 0.2 are 
shown in Figure 5. For each cluster model, the distributions 
tend towards longer t cx when compared to the models using 
the FIR population shown in Figure 4. This is only a conse- 
quence of the larger ft allowed in models that used the FIR 
population, because the bulk of the FIR binary population 
is unobservable. Models that used the FIR population with 
fb = 0.2 had similar or slightly longer f ex than models that 
used the OBS population at the same fb- The difference in 
composition between the FIR and OBS populations, on its 
own, did not impact the amount of time it takes a single BH 
to exchange into a binary. 

Conversely, the composition of the OPT binary popu- 
lation directly impacted i cx and t mrg . Models in cluster D 
that used the OPT population resulted in an average ex- 
change time that was nearly a factor of two smaller than i ex 
in models that used the OBS or FIR populations with the 
same fb- Furthermore, with the OPT population, in 27% of 
the encounters that resulted in a single BH exchanging into 
a binary, the BH exchanged directly into a binary with a 
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Figure 4. Distribution of i GX , the time for a single BH to exchange into a binary. Each panel shows the distribution for models using 
the FIR binary population in cluster B (solid), C (dashed), and D (dotted). Distributions are show for simulations with Mbh = 7 Mq 
and = 0.5 (left), Mbh = 35 Mq and f b = 0.5 (centre), and Mbh = 35 Mq and f b = 0.75 (right). Increasing Mbh does not have a 
large effect on the distribution of t ex . Increasing on the other hand, results in a narrower distribution that peaks at a lower i ox . Note 
that the distributions shown in the centre and right panels include t cx for BHs that became single after merging with a NS, while all 
such systems were promptly ejected from the cluster for the models shown in the left panel. 




10 



logf ex [yr] 



Figure 5. Distribution of t ex for models in cluster B (solid), 
C (dashed), and D (dotted). In these simulations we used the 
OBS binary population, A^bh = 35 Mq, and /(, = 0.2. These 
simulations tend towards much longer t ox despite having a similar 
observable binary population to the simulations shown in the right 
panel of Figure 4. 



NS. Many of these highly eccentric BH-NS binaries quickly 
merged, and such mergers accounted for 10% of all mergers 
in these runs. In models with the other binary populations, 
single BHs exchanging directly into a binary with a NS ac- 
counted for less than 0.3% of all mergers. 

In fact, the only simulation of cluster B that produced a 
run with multiple BH-NS mergers used the OPT population. 



In all other simulations that used cluster B, the post-merger 
BHs did not have enough time to exchange into a new binary, 
gain a NS companion, and merge. Figure 6 illustrates how 
the assumed binary population aflects the merger rate. Each 
model shown used ft = 0.2, cluster D, and Mbh = 35 Mq. 
As was the case with the i cx distributions, the i m r g distri- 
butions for models using the FIR and OBS populations are 
similar. In each, the distribution of initial merger times is bi- 
modal. The first peak corresponds to mergers that occurred 
before the initial BH binary was destroyed in a binary-single 
encounter. The second peak, at t mrg > 3 Gyr is due to BHs 
that lost their companion in a three-body interaction and 
had to re-exchange into a binary system. The model using 
the OPT population does not display this bimodality. The 
OPT population is tuned to interact with the single BH, 
so BHs that lose their companions do not remain single for 
long. In models using all three populations, the distribu- 
tion of subsequent merger times peaks in the largest t mT& 
bin. This is expected because such mergers cannot occur 
until after the BH undergone its first merger. Importantly, 
in addition to increasing the number of mergers that oc- 
curred, and the corresponding merger rate, IZgc, the multi- 
ple merger scenario shifted the peak of all f mrg distributions 
to much later times. This suggests that such mergers could 
be observed during the current epoch. 



As with the single merger scenario, the more concen- 
trated cluster models produced larger merger rates. The i m rg 
distribution for models in clusters B, C, and D are shown 
in Figure 7. Each of the models shown used the FIR bi- 
nary population with /(, = 0.75 and Mbh = 35 Mq. Subse- 
quent mergers dominated the t mrg distribution at t > 1 Gyr 
in cluster D. In the cluster C models, most initial mergers 
did not happen until t>l Gyr, so the subsequent mergers 
were not a significant contribution to the total merger rate 
until late times. In the cluster B model, none of the BHs 
underwent multiple mergers, primarily because many of the 
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Figure 6. Merger time distributions for simulations with = 0.2 using the FIR (left), OPT (centre), and OBS (right) binary population. 
Each panel shows the distribution of t mrg for the initial merger during a run (dashed), subsequent mergers during that run (dotted), 
and the combination of the two (solid). The initial t mrg distributions for the simulations using the FIR and OBS populations show a 
clear bimodality. The peaks above logt mrg = 9.5 in these distributions are due to BHs that lost their binary companions before merging 
with a NS. The bimodality is amplified in the total merger time distributions for these runs by the distinct distributions for initial and 
subsequent mergers. The t mrg distribution in the OPT case is completely dominated by subsequent mergers. 




Figure 7. Distribution of t mrg for simulations using the FIR binary population with = 0.75 in cluster D (left), C (centre), and B 
(right). Each panel shows the distribution of t mrg for the initial merger during a run (dashed), subsequent mergers during that run 
(dotted), and the combination of the two (solid). Subsequent mergers dominate in cluster D. In cluster B, there are no subsequent 
mergers. 



initial mergers did not occur until the BH had evolved in 
the cluster for 6 Gyr. 



Clusters with the highest concentrations and binary 
fractions permitted the largest number of mergers that a 
single BH can take part in. The maximum number of BH- 
NS mergers in a single run, JV max , is listed for each model in 
Table 2. In a run in cluster D using the FIR binary popula- 
tion and ft = 0.75, and in a run in the same cluster model 
using the OPT binary population with ft, — 0.2, a single 
BH merged with 9 NS. These models also resulted in the 
largest proportion of runs that resulted in multiple mergers. 
In each, 73% of the BHs that underwent a single merger, 
underwent at least a second merger. 



4 DISCUSSION 

We have shown that binary-single star interactions in glob- 
ular clusters can produce BH-NS binaries that will merge 
within a Hubble time. The rate at which the mergers oc- 
cur is sensitive to the GC structural parameters, the mass 
of the BH, and properties of the GC's binary population. If 
the BHs in GCs follow a similar mass distribution to those 
found in the disk of the Galaxy, then it is unlikely that 
any high concentration clusters could retain a BH until the 
current epoch. In such clusters, most BHs that managed 
to survive the early period of self-interaction, described in 
previous studies, would gain a NS companion, undergo a 
gravitational radiation driven merger, and be ejected from 
the within ~ 3 Gyr. If, on the other had, GC BHs are more 
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massive than typical field BHs, many clusters could retain 
a BH and that BH could continue to merge with NSs until 
late times. 

4.1 Uncertainties 

The BH-NS merger rates calculated here depend on the un- 
certain distributions of several quantities. The distributions 
of some parameters are better constrained than others (e.g., 
a Abt 1983 and e Duquennoy & Mayor 1991 distributions, 
NS retention fraction Pfahl et al. 2002, and the structure of 
extragalactic GCs Strader et al. 2011). However, both the 
GC BH mass function and the nature of GC binary pop- 
ulations strongly influence the predicted rates and are not 
well constrained by observations. We have sampled several 
plausible scenarios for each, and the results suggest that the 
observed merger rate could be used to probe both the BH 
mass distribution and binary populations present in GCs. 

Some of the simplifying assumptions that we made in 
the models also impacted the predicted merger rates. In our 
calculation of the probability that a single BH would ex- 
change into a new binary, we related the densities of single 
and binary stars using an approximation that diverges for 
large binary fractions, /i,/(l — fb)- In the most extreme case 
considered here, = 0.75, this term increases the encounter 
probability by a factor of 3. Any overestimate due to this 
term is of the same magnitude as errors introduced by other 
approximations made in our calculation. In addition to this 
approximation, in the calculations described above we used 
a relationship between the rates of specific exchange encoun- 
ters to scale the total encounter rate. At the very least, this 
resulted in longer t ex because we only integrated the three- 
body encounter if the probability for a specific exchange 
[(1,2) + BH (1,BH) + 2] was large enough, but we al- 
lowed for any outcome of the three-body interaction. 

We have run two additional simulations to explore how 
these assumptions affected our calculations. In these tests, 
we evolved 500 binaries in cluster D using the FIR binary 
population with /;, = 0.5. In one simulation, marked Retry, 
we scaled the encounter probability using equation 7, as be- 
fore. However, in these runs if the probability of exchange 
was larger than the random number, we re-ran the three- 
body integrator until the desired exchange, (1,2) + BH — > 
(1, BH) + 2, occurred. If the integration resulted in a differ- 
ent outcome (e.g., flyby or an exchange creating a binary 
(BH,2)), we changed the orbital phase of the initial binary 
and the orientation of the encounter, and then simulated 
the interaction again. In another simulation, marked Total, 
we scaled -P(i,bh)+2 by the ratio of the total encounter cross 
sections, not just the single, desired exchange cross section. 
In the mass ratios considered here (~ 1 : 1 : 35), the value 
of ff(X) for the BH to exchange into a binary with the pri- 
mary is roughly equal to a(X) for the BH to exchange into a 
binary with the secondary, and roughly half the value <r(X) 
for a flyby (Sigurdsson & Phinney 1993). This results in a 
factor of 2-3 increase in the encounter probability over the 
probabilities used in the previous runs. For the Total runs, as 
in the runs presented in Section 3.2, if an encounter occurred 
we accepted any outcome of the three-body integration. 

The results of these tests are listed in Table 3. Both 
methods for correcting the probability scaling increased the 
exchange and merger rates. Furthermore, making the correc- 
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Table 3. Probability Scaling Tests 



Class i ox N max TZ-GC 

(yr) (yr -1 ) 

Original 2.01 X 10 9 7 1.7 X 1CT 10 
Retry 9.96 X 10 s 8 2.3 X 10" 10 
Total 1.18 x 10 9 8 2.5 x 10~ 10 

tion either of two ways gave consistent results, illustrating 
that the probability scaling argument used in our models 
is valid. The maximum number of mergers in each of the 
test runs increased from 7 to 8. The fact that this quantity 
did not change by a factor of two, as i ex did, suggests that 
the time it takes a single BH to exchange back into a bi- 
nary may not be the rate limiting step in the BH-NS merger 
process. Finally, these runs suggest that the factor of a few 
uncertainties in the exchange probability resulting from the 
approximations we made can lead to uncertainties of a simi- 
lar magnitude in the BH-NS merger rates. As such, the rates 
given in Table 2 should be treated as having a factor of ~ 2 
uncertainty. 

In addition to uncertainties in the components of the 
models, physical processes and properties of GCs that we 
have omitted from the models could alter the BH-NS merger 
rates. We have not explicitly included an increase in the bi- 
nary density at the cluster centre. Since, next to the BH, 
NS-containing binaries are the heaviest objects in the clus- 
ter they would likely be concentrated near the centre. The 
situation could be similar to the models presented here that 
made use of the OPT binary population, which contained 
a large fraction of NS binaries. These models resulted in 
higher merger rates than models that used OBS or FIR 
binary populations with the same binary fraction. On the 
other hand, a high binary concentration could result in sig- 
nificant heating, expanding the radial profile expected for 
a population of thermalised binaries. A high concentration 
of binaries would also result in binary-binary interactions, 
which were not included in our models. Some binary-binary 
encounters could enhance the merger rate. For example, the 
interaction of a BH-MS binary with a NS-MS binary could 
produce a BH-NS binary. However, there is a high likelihood 
that a binary-binary encounter will result in the disruption 
of a binary. Because the rate of binary-binary encounters 
increases with fb, this limits the binary fraction. So, binary- 
binary interactions could cause a reduction in the BH-NS 
merger rate, compared to the rates predicted here, by limit- 
ing fb and making it more difficult for a single BH to find a 
new binary companion. Detailed models that include binary- 
binary interactions are necessary to determine which process 
dominates. 

4.2 LIGO detection rate 

We estimate the LIGO detection rate following 
Belczynski et al. (2007) and Banerjee et al. (2010). We 
computed the detection rate with 

Hligo = %Pgc(dA mch ) 5/b ) TZgc, (8) 

where pGC is the space density of GCs, m c h is the chirp 
mass of the binary (m c h = (mim2) 3//5 /(mi + ma) 1 ' 5 ), and 
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m c h,NSNS is the chirp mass of a binary with two 1.4 M© 
NSs. Do is the maximum distance from which the gravita- 
tional waves from a merging NS-NS binary can be detected, 
and has a value of 18 and 300 Mpc for LIGO and Advanced 
LIGO (AdLIGO), respectively. Using the galaxy luminosity 
function parameters estimated in Croton et al. (2005) and 
the specific frequency of GCs per 8.5 x 10 7 L© of galaxy 
luminosity presented in Brodie & Strader (2006), we calcu- 
lated pgc = 8 Mpc -3 . Our calculations above showed that 
IZgc is strongly dependent on the cluster mass. To compute 
a merger rate that is averaged over all GCs, IZgc, we used 
the GC mass function presented in McLaughlin & Pudritz 
(1996). 



Although we have shown that the bulk of such mergers 
occur quickly, and therefore at high redshift, we compute the 
LIGO detection rate for models with 7 M© BHs for illustra- 
tive purposes. For these BH-NS binaries, m c h = 2.6 Mq. 
For the models that do not allow single BHs, we found 
TZgc = 3.4 x 10 -12 yr -1 GC -1 . To make a conservative 
estimate for the AdLIGO detection rate, we assumed that 
10% of clusters retain one BH and arrived at the rate of 
2 x 10~ 3 yr -1 . Using instead the rates from runs that did 
allow single BHs to exchange back into binaries, we found 
an AdLIGO detection rate of 10~ 2 yr -1 . Even if mergers 
between 7 M© BHs and NS were occurring in GCs at the 
current epoch, the event rates are so low that they would not 
be detected by AdLIGO. For example, the rates calculated 
for cluster models B and C may be applicable to massive in- 
termediate age clusters. With ages of 1-5 Gyr, these clusters 
are young enough that any retained BHs would only now be 
merging with a NS. However, the space density of massive 
intermediate age clusters clusters is even lower than that of 
GCs, so the conclusion remains the same. If GC BHs fol- 
low the same mass distribution as field BHs, BH-NS merg- 
ers resulting from binary-single star interactions in clusters 
cannot be detected by AdLIGO. 



4.3 Comparison with previous studies 

Previous studies have estimated the BH-NS merger rate 
in the field and in GCs. The field rates are typically cal- 
culated using binary population synthesis models that are 
subject to uncertainties in common envelope evolution, 
wind mass loss from high mass stars, and natal kicks, re- 
sulting in a wide range of predicted AdLIGO detection 
rates, 0.68-42.8 yr _1 (Sipior & Sigurdsson 2002; Pfahl et al. 
2005; Belczynski et al. 2007; O'Shaughnessy et al. 2010; 
Belczynski et al. 2010). However, a case study of the likely 
BH-NS progenitor Cyg X-l presented in Belczynski et al. 
(2011) avoids many of the poorly constrained components 
of binary population synthesis and predicts much lower field 
BH-NS merger detection rates in the range (0.4 — 2.8) x 
lO" 2 yr- 1 . 

Our results are consistent with previous studies of 
the GC merger rate that conclude the cluster merger 
rate for both NS-NS and BH-NS binaries is much lower 
than the field rate for corresponding pessimistic, realis- 
tic, or optimistic predictions (Phinney 1991; Grindlay et al. 
2006; Sadowski et al. 2008). However, if the uncertainties 
in the field merger rate push it towards the pessimistic 
estimates, which is suggested by the models presented in 
Belczynski et al. (2011) that circumvent many of these un- 
certainties, and multiple mergers increase the GC rate, 
then the BH-NS mergers detected by AdLIGO might be 
dominated by, or comprised entirely of, systems that form 
through multibody interactions in GCs. 
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The detection rate could be much higher than this 
conservative estimate if GC BHs are more massive than 
typical field BHs. Higher mass BHs enhance the detec- 
tion rate in two ways. First, more massive binaries can be 
detected out to larger distances. Second, since the post- 
merger kick is suppressed for high mass ratio inspirals, 
the BH remains in the cluster and can undergo multiple 
mergers with NSs. Models with the FIR binary population 
and ft = 0.75 produced the largest averaged merger rate, 
TZgc = 4.2 x 10~ n yr" 1 GC _1 . Assuming that 10% of GCs 
retain a 35 M© BH, this corresponds to a AdLIGO detec- 
tion rate of 0.14 yr -1 . If 50% of GC retain a BH the rate 
would approach 0.7 yr -1 . The FIR binary population with 
ft = 0.75 represents an optimistic, yet physically motivated 
population. On the low merger rate side, models with the 
OBS population and fb = 0.05 would yield an AdLIGO 
detection rate of 0.04 yr -1 if 10% of clusters retained a 
35 M© BH. Although the AdLIGO detection rates are not 
promising, the predicted merger rates are large enough that 
BH-NS mergers in GCs will be seen in the next generation 
of gravitational wave detectors (e.g., the Einstein Telescope 
Punturo et al. 2010). 
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